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We investigate the properties of conduction electrons in single-walled armchair carbon nanotubes 
(SWNT) in the presence of both transverse electric and magnetic fields. We find that these fields 
provide a controlled means of tuning low-energy band structure properties such as inducing gaps 
in the spectrum, breaking various symmetries and altering the Fermi velocities. We show that 
the fields can strongly affect electron-electron interaction, yielding tunable Luttinger liquid physics, 
the possibility of spin-charge-band separation, and a competition between spin-density-wave and 
charge-density-wave order. For short tubes, the fields can alter boundary conditions and associated 
single-particle level spacings as well as quantum dot behavior. 



I. INTRODUCTION 

The astounding range of experimental and theoretical 
studies performed on carbon nanotubes [l[ has revealed 
a spectrum of physics characteristic of strongly corre- 
lated low-dimensional electronic systems [ll . The under- 
lying graphene lattice structure of these tubes uniquely 
affects band structure, effective dimensionality, Coulomb 
interaction effects, and the quantum dot behavior ex- 
hibited by short nanotube segments. The band struc- 
ture shows differing behavior depending on various fac- 
tors such as chirality, applied gate potentials, boundary 
conditions at the tube ends and mechanical stress [1, 
In single- walled armchair nanotubes (SWNT) which 
are the entities of interest here, gapless linearly dispers- 
ing modes endow the nanotube with its peculiar quan- 
tum wire properties. As described theoretically and as- 
certained experimentally, interactions within the modes 
of this effectively one-dimensional system cause it to be- 
have as a Luttinger liquid characterized by non-Ohmic 
conductances 0, [g, 0, II, 9 • Tubes placed between tun- 
nel barriers act as quantum dots (lol . [ll|, which, while 
displaying zero dimensional physics such as Coulomb 
blockade behavior, retain some higher dimensional traits 
such as hosting plasmons typical of one dimension and 
band degrees of freedom attributed to the underlying 
graphene lattice. Potentially invaluable to applications, 
these nanotube quantum dots have been proposed as el- 
ements of quantum devices and the quantum states of 
blockaded electrons have been regarded as candidates for 
units of quantum information |12l | . In each of these as- 
pects, the presence of applied fields can dramatically al- 
ter the nanotube's behavior; here we present an extensive 
study of the effects of electric and magnetic fields applied 
transversally to the axis of the nanotube. 

At the level of the band structure, it has been shown 
that a parallel magnetic field can have the striking effect 
of converting a metallic tube to a semiconducting one 
by way of inducing a gap [l^l, and vice- versa, an effect 
discernible in conductance, Coulomb blockade and scan- 



ning tunneling microscope (STM) measurements. Here, 
instead of a parallel field, we discuss transverse field con- 
figurations (both electric and magnetic) and the condi- 
tions under which a band gap opens up or the spectrum 
remains gapless in armchair SWNTs. In the latter case, 
we demonstrate, via band-structure calculations, simul- 
taneous breaking of the valley degeneracy (of the two dis- 
tinct Dirac points), the left-right-mover degeneracy, and 
the particle-hole symmetry. Moreover, the fields yield a 
non-negligible reduction in the Fermi velocity of conduc- 
tion electrons traveling along the tube. We show that for 
certain configurations of fields, the ground state of the 
tube can even be made to carry finite current. 

Transverse fields provide an excellent means of alter- 
ing the ratio of interaction strength to the Fermi energy 
in SWNTs. This makes nanotubes potentially the only 
systems to date in which the associated Luttinger liquid 
physics can be tuned in a controlled fashion. As described 
in previous work, either an electric field or a magnetic 
field 

alone suffices to change the value of the Lut- 
tinger liquid parameter from that measured in field-free 
environments. The magnitude of the electric fields re- 
quired to bring about a significant change are well within 
current experimental reach [l3|. Here, we find that our 
approach reproduces these results. We show that such 
a tuning of Luttinger parameters can mediate a transi- 
tion from the system showing tendencies towards spin- 
density-wave (SDW) ordering to that of charge-density- 
wave (CDW) ordering. Furthermore, in addition to the 
tuning of the Luttinger parameter presented in previous 
works for the net charge density [ij, E^j [3 1 Siid that 
Luttinger-type interactions become manifest in modes as- 
sociated with the density differences between nanotube 
bands as well. Thus we predict that akin to spin-charge 
separation, transverse fields can induce a spin-charge- 
band separation wherein the three degrees of freedom 
move at different velocities. 

The above results are discussed in the case of an in- 
finite system. For short tubes or finite length segments 
formed by tunnel barriers, boundary effects need to be 
taken into account. We find that applied fields influ- 




FIG. 1: A (5,5) carbon nanotube in the presence of transverse 
magnetic (pointing in the —y) and electric fields. The carbon 
atoms belonging to the A and B sublattices are indicated by 
dark (blue) and light (green) shading, respectively. 

encc multiple aspects of sliort nanotubes. First, fields 
can alter the single-particle energy level spacing of the 
tube. Here, we carefully account for the effect of the tube 
ends in the case that the left and right movers travel at 
different speeds. Second, the charging energies become 
field dependent, and third, the plasmon spectrum varies 
in accordance with the first two effects. We have con- 
ducted a comprehensive analysis of the short nanotube 
as a finite-sized Luttinger liquid and show how all three 
effects can be captured. Within this description, wc dis- 
cuss the structure of field-dependent Coulomb blockade 
peaks and how the presence of both electric and magnetic 
fields acts as a means of manipulating quantum states of 
the dot (the effect of a magnetic field alone has recently 
been discussed by Bellucci and Onorato [isj). 

The outline of this paper is as follows. In section II wc 
present the formulation and results of our band struc- 
ture calculation. In section III we formulate an effective 
one-dimensional Hamiltonian which takes into account 
field effects. We bosonize this Hamiltonian, describing 
interaction effects in terms of Luttinger liquid physics. 
In section IV we investigate the various Luttinger liquid 
phases and the feasibility of using fields to access such 
phases. In section V we discuss field-tuned quantum dot 
physics. Finally, in section VI we present the highlights 
of our results and discuss their relevance to experiments. 

II. BAND STRUCTURE IN TRANSVERSE 
FIELDS 

We briefiy recapitulate the band structure of an in- 
finitely long armchair tube in the absence of any fields Q ■ 
The electronic properties of graphene are well described 
by a tight-binding model in which electrons hop be- 
tween nearest neighbors of the underlying hexagonal bi- 
partite lattice (sublattices here labeled A and B) with 



an associated energy (the hopping integral) of < « 3 
eV. An armchair carbon nanotube can be regarded as 
a sheet of graphene rolled along the (n, n) -direction (for 
notation of the chirality, see e.g. Ref. Q), denoted by 
s. This gives rise to states of quantized momentum 
ks = (0, 27r/L, 27r(2r? — 1)/L) where L = \/3na is the 
circumference of the tube, and a = VSac where Uc ~ 0.15 
nm is the nearest carbon-carbon distance. The result- 
ing series of one dimensional bands can be described by 
the wavevectors k = {k^ , ^ ) where kx is the quasimo- 
mentum parallel to the tube's axis and M is the state's 
angular momentum about the tube's circumference. A 
convenient set of basis states is given by the following 
linear combination of atomic orbitals 

R<£A/B 

where \R) is the 7r-electronic state of the atom located at 
R. The sum runs over the n atoms in the unit cell that 
belong to either the A or B sublattice. At half filling, 
the associated dispersion has low energy excitations near 
the so-called Dirac points of the form e = zthvp\k — akpl 

where vf = « 8 x 10^ m/s, kp = 47r/3a and a = ±. 
Thus, k = akp label the two inequivalent Fermi points. 

The setup of interest is shown in Fig. [1] An external 
magnetic field is applied in the negative y-direction; an 
applied (transverse) electric field makes an angle x with 
the magnetic field. These fields give rise to scalar and 
vector potentials 

U{s) = |e|i?cos(-i-x) (2) 

A = -Bzx, (3) 

respectively (i? ~ L/2tt). The x-axis runs parallel to the 
tube's axis and the additional coordinate s measures the 
circumferential distance starting from the negative y-axis 
(a positive value of s corresponds to a counterclockwise 
rotation as one looks along the x-axis in the positive di- 
rection). 

These external potentials are easily accommodated 
within the tight binding approach. In Eq. ^ wc have 
selected a gauge that is independent of x and thus kx 
remains a good quantum number. The hopping matrix 
elements in the presence of the fields are given by 

\H\ = - g«(fe-fl-fc'--R')+Tf (Gh-G'«) 

ReB.R'eA 

where the sum runs over nearest neighbors i?, R' and 

G, - « / dX in - rj) ■ A{t'+ A(f, - rjj) (5) 
Jo 

is the Aharonov-Bohm phase associated with the mag- 
netic field Q . The dimensionless parameter b is given by 

b = B ^l^}^ ■ Numerically, for an {n,n) nanotube, the 
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FIG. 2: Spectrum of a (5, 5) carbon nanotube near the a = + 
Dirac Fermi point (field-free value kp = a4n/3a indicated by 
the vertical line) in the presence of an external perpendicular 
field with Uo/t = (black dotted), 0.2 (red dot-dashed), 0.4 
(blue solid) as the crossing moves to the right. The horizontal 
axis indicates the value of q, where q = ka; the vertical axis 
is given in units of t, the hopping integral (t fa 3 eV). 



FIG. 3: Spectrum of a (5, 5) carbon nanotube near the a — + 
Dirac Fermi point (field-free value kp = a4n/3a indicated by 
the vertical line) in the presence of a magnetic field b = 
(black dotted), 0.2 (red dot-dashed), 0.4 (blue solid) (as the 
crossings move to the right). The horizontal axis indicates 
the value of q, where q — ka; the vertical axis is given in units 
of t, the hopping integral (t ~ 3 eV). 



magnetic field in Teslas is related to the dimensionless 
parameter b via B « 8.1 x 10^ x b/n'^ T. The scalar po- 
tential gives rise to an on-site potential described by the 
matrix element 

($1 \H\ $i) = \H\ <fi) = ^e^^^ (6) 

for e' = £ ^ I mod n where U ~ \e\ER/t. The electric 
field strength in V/nm is related to U by E 42 U/n 
V/nm for a tube with chiral vector {n,n). That these 
matrix elements mix states of different angular momen- 
tum has a straightforward classical analog: a charged 
particle on the surface of a cylinder will circulate around 
its circumference as a result of the applied fields 

We have studied the effects of the fields perturbatively 
in b and U. In the vicinity of the Fermi points, the 
left and right moving bands are nearly degenerate, so 
care must be taken in applying perturbation theory. The 
details of this calculation are presented in Appendix |^ 
where we carry out perturbation theory to second or- 
der. There are three cases of interest that we summarize 
below. An illustration of these three cases is shown in 
Figs, m [21 m andO While some of the effects discussed 
are rather small for standard SWNTs, we note that our 
band structure analysis can be applied to multi-walled 
nanotubes as well in which case a larger radius yields 
more pronoimced effects. 



1. Case ofE = or B = 

For a single field, the most salient features of our band 
structure calculation are the reduction in Fermi velocity 
and the shift in Fermi momentum. Semiclassically, the 
reduction in Fermi velocity can be ascribed to the deflec- 
tion of the electrons by the fields leading to a reduction 



in the velocity component along the tube. Furthermore 
the bands remain gapless. 

For a magnetic field, we find a reduced Fermi velocity 
given by 

VF =vf{1- Avib^) . (7) 

The first order correction in b vanishes because vp must 
be an even function of b by symmetry. The term Avi 
is a function of n and is given by Eq. (jA4bp . This term 
depends on the geometric details of an armchair tube; for 
large tubes Avi « 1/3. An experimental observation of 
such a reduction will require strong fields. For example, 
a (20,20) armchair tube with a 20 T field gives a 0.4% 
reduction in the Fermi velocity. 
For an electric field we have 

VF = WF (1 - Awsf/^) . (8) 

A field of strength 0.1 V/nm corresponds to a reduction 
in the Fermi velocity of a (10,10) tube of roughly 10%. 
The term Av2 ~ r? jiv^ for large tubes; its exact form is 
given by Eq. \k4c} . 

As mentioned in the introduction, the reduction of the 
Fermi velocity has been noted by several authors [3, [H, 
[iGj . In these papers the carbon nanotube was modeled 
as a smooth cylinder and the low-energy electronic be- 
havior was put in by hand. Our results are in agreement 
with these results in the limit of small fields and large 
n (in the regime that perturbation is valid). Addition- 
ally, by taking into account the geometry of the armchair 
nanotube we find that either a magnetic or electric field 
alone will shift the Fermi points. That is, the nanotube 
still has the same low energy spectrum with renormalized 
values oikp (whose precise form is given by Eq. (|A5[) ). 
The band structure of a (5,5) nanotube in the presence of 
electric and magnetic fields of various strength is shown 
in Figs. [2] and [21 respectively. 
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FIG. 4: Spectrum of a (5, 5) carbon nanotube near a — ~ 
Dirac Fermi point (field- free value akp = — 47r/3a indicated 
by the vertical dashed line) in the presence of transverse elec- 
tric and magnetic field (U/t = 0.2 and b — 0.4). The angle be- 
tween E and B being (black dotted), 7r/4 (red dot-dashed), 
7r/2 (blue solid) (from outer to inner). The horizontal axis 
indicates the value of q, where q = ka; the vertical axis is 
given in units of t, the hopping integral (t « 3 eV). 
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to expect the band structure will mimic the behavior 
of a classical velocity selector. Indeed the asymmetry 
of the velocity in the right- and left-moving branches is 
expected from a elementary consideration. For a charge 
particle moving on a smooth cylinder in the presence of 
mutually perpendicular transverse magnetic and electric 
fields, the direction of the force caused by the magnetic 
field depends on whether the particle is moving along 
the tube in one direction or another, whereas the electric 
field remains the same. Hence, for fixed kinetic energy, 
the force in the transverse direction causing the particle 
to spiral is different for different directions of motion, and 
giving rise to different velocities along the axial direction. 

Another prominent feature of the band structure in 
this case is a relative energy shift with respect to the two 
Fermi points. Near the two Dirac points kp 
have 



±ft we 

3a 



t{k) = hrv. 



atAs + 0{Ak^), (10) 



where a (a = ± for fcp w ±|^). By definition, kp is the 
momentum for which the left and right moving bands for 
a given Fermi point are degenerate and again is generally 
different from 47r/3a in the presence of fields. The precise 
forms oi kp and As are given by Eqs. (|A5P and (jA6[) . 
respectively. The solid blue lines in Figs. [H (for q = -|-) 
andO (a = — ) indicate low-energy dispersion of a (5, 5) 
nanotube in the presence of crossed electric {v = 0.2) and 
magnetic (b = 0.4) fields. 
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FIG. 5: Spectrum of a (5, 5) carbon nanotube near a = -I- 
Dirac Fermi point (field- free value akp = 47r/3a by the ver- 
tical line) in the presence of transverse electric and magnetic 
field [U/t = 0.2 and b = 0.4). The angle between E and B 
being (black dotted), 7r/4 (red dot-dashed), 7r/2 (blue solid) 
(from outer to inner). The horizontal axis indicates the value 
of q, where q = ka; the vertical axis is given in units of t, the 
hopping integral (t « 3 eV). 



2. Case of ELB (x = t^/2) 

Mutually perpendicular fields break both the time- 
reversal and particle-hole symmetry of the band struc- 
ture. The left and right movers now move with different 
speeds. For a magnetic field in the negative y-direction 
and an electric field in the positive z-direction we have 



Vp 



(9) 



where r is -I- for right movers and — for left movers. The 
expressions for Av^ is given by Eq. (jA4d[) : for large tubes 
Alia « tz/tt. For the fields we consider, E/B is roughly 
the same order of magnitude as vp and thus it is natural 



The presence of both electric and magnetic fields will 
generically open gaps at the Fermi points. The size of 
this gap is given by 



2V3( 
i tbU—— 
1 + 2 



' o cos , 

^ COS Y 

-2cos:f ' 

3n 



(11) 



For example, a (15,15) tube parallel electric and magnetic 
fields 1 V/nm and 10 T respectively gives egap/ks ~ 
4 meV. For electric and magnetic fields which are not 
parallel, this gap is weakly indirect. 

The degeneracy associated with the graphene Fermi 
points stems from the equivalence of the two sublattices. 
For example, a two-dimensional honeycomb lattice with 
A and B sublattices composed of different types of atoms 
would generically have a gap Now, in a nanotube, 
gaps can arise for different reasons. For example, the 
gap associated with a semiconducting tube occurs be- 
cause the quantized bands miss the Dirac point of the 
underlying graphene lattice. Such a gap can be closed by 
applying a magnetic field (of a specific strength) parallel 
to the tube's axis. In the present case however, the gap 
arises from an energy difference associated with the two 
sublattices and therefore can not be so closed. Not sur- 
prisingly, the gap described here vanishes precisely when 
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the electric and magnetic forces on a classical charged 
particle traveling along the tube are either parallel or 
antiparallel. 

From Figs. S] and [5] we clearly see that as the angle 
between the fields varies, the particle-hole symmetry and 
the valley degeneracy, as well as the symmetry between 
left- and right-moving spectra, are broken. The degree of 
breaking of these symmetries is greatest in the case that 
the fields are perpendicular. On the other hand, these 
symmetries are preserved when the fields are parallel or 
anti-parallel and the gap is at its maximum. This angle- 
dependent gap can be manifest in the transport mea- 
surement, e.g., the conductance can be 4e^//i, or 2e^/h 
or zero, depending on whether there are 4, 2 or zero con- 
ducting channels (including spin). 

4- Low-energy modes 

In order to study the low-energy physics of electronic 
excitations about the Fermi energy for all the cases de- 
scribed above, we can identify the fermionic operator on 
the cylindrical surface of the tube as approximately 

*cr(a:, s) = ^ ifpaix, S)'lppaa{x), (12) 
pa 

where tppaaix) is the one-dimensional field operator at 
the point x along the tube axis associated with the 
A and B sublattices {p ~ A{+)/ B{ — )), Fermi points 
kp ~ aAir/Sa (a — ±), and spin a (-1-)/ | (— ). 
The Bloch wavefunction Lppa{x, s) retains detailed infor- 
mation about the response of the electronic wave func- 
tions to the applied fields [1, 0]. While the sublattice 
basis (indexed by p) is convenient for discussing these 
wavef unctions, it does not diagonalize the hopping in- 
teraction. We therefore transform to the right and left 
moving basis (r — R(+) / L(—)) via the transformation 
V'pcKT = Upripracr whcrc U^ayU ^ [31. This gives a 
kinetic energy term 

Ho = -ih'^ J dx rvri^l^d^^jra, (13) 

where, as shown above a generic field configuration can 
give rise to the possibility that vr ^ v^. 

The gap which arises as a result of the presence of both 
electric and magnetic fields [Eq. pTjl ] can be incorporated 
into our Hamiltonian by a mass term 

i?gapl - E / (14) 
racr 

Note that this mass term is of a different form than that 
considered in the exhaustive zero-field study by Egger et 
al.Q (which had the form i^^^ripl^^tl^-racr)- Similarly, 
the effect of the shift associated with the Fermi points 
can be described by the term 

-ffgap2 [ dx a^VL^V-rao-. (15) 



III. LUTTINGER LIQUID FORMULATION 

In this section we specialize to the case of gapless 
modes for which interactions can be easily incorporated. 
As the simplest case, when only a magnetic or electric 
field is present, the kinetic piece of the Hamiltonian given 
in Eq. (|13p now has ~ = vp and egap = 0, i.e.. 

Ho = -ihvp'Y^ j dx ipl^^d^ipraa- (16) 

racr 

We use standard approaches such as bosonization to 
study the effect of fields. Where appropriate, we include 
discussions for the asymmetric case of ^ v^. We 
closely follow the approach of Ref. [3] whose lucid peda- 
gogical exposition we do not repeat but instead confine 
our discussion to the new field-dependant features. 

The presence of fields does not alter the fact that the 
electrons on a tube are locked into their lowest energy 
radial mode. Hence, it is possible to study the low en- 
ergy excitations in the presence of interactions using an 
effective bosonized ID Hamiltonian. Bosonization offers 
a great simplification since many of the quartic interac- 
tion terms in the fermionic language become quadratic 
once they are bosonized. 

The interaction term takes the general form 

H,r^t = ^JdrJdv' *t(i.)^t,(r')C/(r-r')*.'(r')*^(r), 

(17) 

where ^„ (r) is the field in Eq. p2p describing low-energy 
electrons. 

Following 0, we employ the form of the Coulomb in- 
teraction on the surface of a cylinder given by 

U{x -x',s- s') = ^^^'^ =, 

^{x - x')2 + 4i?2 sin2 (£^) + al 

(18) 

where R is the radius of the tube and x and s the co- 
ordinates defined in the previous section and Uz ~ a is 
roughly the thickness of the graphene sheet . The form 
of the interaction in Eq. (|17|) is explicitly given in terms of 
two-dimensional integrals. An effective interaction term 
involving purely linear integrals along the tube's axis can 
be obtained by expressing 5'(j(r) in terms of the linear 
and circumferential fields as in Eq. (|12p and integrating 
out the circumferential degrees of freedom from Eq. (|17p . 
Here we remark that integration over the circumferen- 
tial degrees of freedom generates field-induced terms that 
give rise to novel physics and can be traced to the depen- 
dence of the ground state wave function on non-zero an- 
gular momentum states in the circumferential direction. 

The resulting effective interaction involves two-particle 
scattering processes between fermions moving along the 
tube axis denoted by the fields ipraa- The associated scat- 
tering processes can be classified according to whether 
the incoming particles preserve their Fermi point quan- 
tum number a when scattered, - forward scattering 
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{aFSy, or scatter across the Fermi surface - backscatter- 
ing (aBS). In Ref. 0, a further distinction is made in 
the forward scattering events depending on whether the 
interaction potential is homogeneous over the circumfer- 
ence of the tube (aFSO) or is able to distinguish the 
microscopic differences which arise between sublattices 
{aFSl). At this point, in order to obtain an effective 
low-energy description of the interacting electrons, the 
one-dimensional fermionic operators can be bosonized as 
in Ref. 0]: 



: exp [iakpX + i^Praa] ■ 



(19) 



where 



+ raOc- (20) 



— (0C+ +r6'c+ + a 
acj). 



The bosonic fields Lp^s satisfy the commutation relations 

VPraa{x),(Pr'a'cj'{x')] = -i7rr(5rr"5^^' sgn(a; - x') , (21) 

where r = ± denotes the left- and right-movers, a = ± 
indicates the Fermi points, and cr = ± represents the spin 
direction (t / i). The the so-called Klein fac- 

tors; they enforce the anticommutation relations between 
different channels, 



(22) 



Moreover, the fields 9js{xys (with j = c/s and S = ±.) 
and their dual fields (jfjs (x) [both are linear combinations 
of ifrairix)] in turn satisfy 

[(l)js{x),9.j>s'{x')] = --Sjj'Ss5>sgn{x - x') (23) 
The effective density in a given channel takes the form 



Praa{x^ — dx^racri^X^. 

Ztt 



(24) 



The kinetic energy associated with the linearly dispers- 
ing fermionic modes is quadratic in the bosonized fields. 
As for the interactions, the dominant contributions also 
come from quadratic terms reflecting net density-density 
type interactions. As in the field-free case, the aFSO pro- 
cess has one such contribution from the usual Coulomb 
form involving the net charge density, which in the 
bosonized representation is given by 

Hc.FS,o = ^ fdxV{k^ 0) {d.J,+f , (25) 



where 



V{k) 



2e' 



(|lnfci?| +co) 



(26) 



is the Fourier transform of V{x) and cq is a function of 
n (see Eq. (|B2l) ). 



The presence of either an electric or magnetic field 
gives rise to additional quadratic terms in the aFSO pro- 
cess . These terms have their origin in the non-zero angu- 
lar momentum components of the circumferential wave- 
function Lpraix,s) in Eq. (jl2[) . A detailed accounting of 
the radial wave functions (see appendix B) shows that an 
electric field contributes a term 



H. 



aFS,E 



dx 



(27) 



whereas a magnetic field provides a contribution 

H^FS,B = J dx (^^^ (^^^ h^b^ {dxe.^f . (28) 

In the presence of mutually perpendicular transverse elec- 
tric and magnetic fields, there is an additional contribu- 
tion to the interaction 

HaFS,BE = (j^^ h'iUh{dxOc-){dx(t)c+)- 

(29) 

The values of /ii , /12, and /13 are given in appendix B. The 
Luttinger liquid Hamiltonian takes into account all these 
quadratic terms; other terms emerging from the interac- 
tion potential are sub-dominant and can be considered 
perturbatively. 

An intuitive physical picture for the origin of the field- 
dependant interaction terms above can be obtained by 
noting that in Eq. (|A7p the charge density for p^a is of 
the form 1 + tiru + t2ab. This is broadly consistent with 
the magnetic field coupling to momentum via (p — eA/c)^ 
and hence to akp, resulting in a term ~ ah. Moreover, 
the electric potential differs slightly between adjacent A 
and B sites, and in turn couples to left- and right-movers 
differently, resulting in a term ^ rU . Equation (fT7|) 
thus gives rise to interactions among the charge densi- 
ties dxOc± and current d(j)c+- We therefore obtain the 
following nonvanishing terms: (1) dxdc+, (2) dxOc- ~ b, 
and (3) dx4>c+ ~ u. Hence, one expects the resulting 
interaction terms given above. 

The sum of the kinetic energy and interactions de- 
scribed by Eqs. (HI]), (gS]), and ^ gives the total 
Hamiltonian that we focus on in this paper 



Htot — Hq + HaFS.O + HaFS.E + HaFS,B- 



(30) 



For the case of a single field (either magnetic or electric) 
the dispersion remains symmetric and therefore Eq. (j30p 
can be written in the form 



EVa 
2 

a—±c/ s 



dx 



Fa 



Ka idx4>aY 



(31) 

Ka = 1 reflects no interactions in the 'a' sector and Ka < 
1 reflects repulsive interactions. For the asymmetric case, 
the diagonalization of the Luttinger liquid Hamiltonian 
is technically more complicated but conceptually simple 
given the quadratic form of the relevant terms [l9| . Below 
we discuss the form of the various Luttinger parameters 
and the related physics. 
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FIG. 6; Luttinger parameter Kc+ of a 500 nm long nanotube 
as a function of electric field strength E (in V/nm) for dif- 
ferent values of n: n = 10 (solid blue line), n = 15 (red 
dot-dashed line), and n = 20 (black dotted line). 



FIG. 7: Luttinger parameter Kc- of a 500 nm long nanotube 
as a function of magnetic field strength B (in Teslas) for dif- 
ferent values of n: n = 10 (solid blue line), n = 15 (red 
dot-dashed line), and n = 20 (black dotted line). 



A. Luttinger parameters for B = Q, E ^ 

As seen above, for just an electric field present, interac- 
tions only affect the net charge sector c+ and rcnormahzc 
the velocity in this sector 



8e 



(32) 



Vc+/Kc+ = VF + (|lnfc,i?| + co), (33) 

TTK 

where kc ~ 1/ L is the lower cutoff provided by the length 
L of the tube. This yields a Luttinger parameter value 



Kc+ = 



\ 



■kk,vf{U) 



■kkvf (U) 



|lnfc,i?| +Co) 



(34) 



{Kc+ + Kc+ ~ 2)/8 for tunneling into the end or 



in the low-field limit, where vf{U) is the field reduced 
Fermi velocity calculated in the previous section. Fig- 
ure [6] shows the dependence of Kc+ on U. In the 
other sectors, the Luttinger parameters retain their non- 
interacting value Kc- = Ks± = 1- The Luttinger model 
predicts power law behavior for the tunneling density 
of states with exponents Uend = iKc+ ^ l)/4 and 

Oibulk 

bulk of a tube respectively. 

A smaller value of Kc+ implies stronger repulsive in- 
teractions in the net charge sector. Here, two tendencies 
compete; in Eq. (p4|) . vf{U) < vf increases the inter- 
action strength relative to the kinetic energy. On the 
other hand, the direct effect of the field, as reflected in 
the hiU'^ term in Eq. (|34p is to decrease the relative in- 
teraction strength. The former effect is dominant of the 
latter, hence as shown in Fig. [51 Kc+ is monotonically 
decreasing in increasing field strength. 



B. Luttinger parameters for B ^ 0, E = 

In the presence of a magnetic field, the Luttinger pa- 
rameter deviates from the non-interacting value of unity 
not only in the net charge sector. 




nK,VF{b) 

but also in the relative charge sector 



{\\iikcR\ + co), (35) 




TTKVFib) 



(36) 



The spin sectors remain unaffected; Ks± = 1. 

For both cases c+ and c— , the reduction in the Fermi 
velocity vpib) < vf increases the relative interaction 
strength, thereby decreasing the values of Kc±. In ad- 
dition, as reflected in the /i2&^ coefficient in Eq. (|36p . 
the field directly decreases Kc~ via the interaction term 
of Eq. (pS)) discussed above. Figure [7] shows the depen- 
dence of Kc- on the magnetic field B. In this case, the 
Luttinger model predicts power law behavior for the tun- 
neling density of states with exponents agnd — {Kc+ + 
K~_} - 2)/4 and abuik = {Kc+ + K'^ + Kc- + K'} - 4)/8. 
Unfortunately, this effect is small; even for fields as large 
as lOOT, Kc- is reduced from unity only by about one 
half of a percent for a (20,20) tube (see Fig. [T]). Nev- 
ertheless, for multi-wall nanotubcs, the outermost layer 
can have n as large as 100, and hence, the field only needs 
to be as large as 4T to see a 0.5 % reduction. 

The deviation of Kc- from unity leads to the fasci- 
nating prospect of spin-band-charge separation. In one- 
dimensional systems, the possibility of spin-charge sep- 
aration stemming from different interactions within the 
two sectors has been extensively discussed and observed 
in the case of etched quantum wires [20| . Here we predict 
that nanotubes in transverse magnetic fields can undergo 
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FIG. 8: Luttinger parameter Kc+ of an n = 10 (solid blue 
line), n = 15 (red dot-dashed line), and n — 20 (black dotted 
line) carbon nanotube as a function of electric field strength 
E (in V/nm) in the presence of a 5T magnetic field (tube has 
K°+ = 0.2). See Eq. §SI 
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FIG. 9: Luttinger parameter Kc+ of an n = 10 (solid blue 
line), n = 15 (red dot-dashed line), and n — 20 (black dotted 
line) carbon nanotube as a function of electric field strength 
E (in V/nm) in the presence of a 5T magnetic field (tube has 
= 0.2). See Eq. §Si 



yet another separation due to the interactions in the c— 
sector. Thus we propose that in this case, the four modes 
travel at three different velocities, Vc+ = vp{b)/ Kc+, 
Vc- = VF{b)/Kc- and Vs± = VF{b). 



C. Luttinger parameters for EA-B 

As shown in Fig.[7l the spin-band separation discussed 
above for purely electric fields is a small effect even for 
very large fields . The prospect of observing this effect is 
greatly improved for crossed electric and magnetic fields. 
In this case, the term which mixes the c+ and c— sectors 
is given by Eq. (^5]) and the values of /13 are significantly 
larger than the corresponding terms hi and /i2 (see ap- 
pendix B). As mentioned previously, the case of crossed 
fields gives rise to an asymmetric dispersion, and it leads 
to cross terms {dxOa){dx4'a) with coefficients proportional 



to the difference vn — v^, which is quadratic in the ex- 
ternal fields. The full treatment of these and the cross 
term in Eq. (^5]) is beyond the scope of the current paper. 
However, it is possible to get a sense of the order of mag- 
nitudes involved by noting that in the present case the 
dominant field effect is given by Eq. We shall in the 
following study the effect of this cross term by ignoring 
the asymmetry in the dispersion. The relevant terms in 
the Hamiltonian are 



H = 



VF 



dx 



dx 



1 



+9 



J dx{dxOc-){da;(f>c+), 



(37) 



with 



vp/v^^ defining the field- free values. In di- 
agonalizing the resulting Hamiltonian, care is required to 
ensure that transformed fields respect commutation re- 
lations such as those of Eq. pT|) . It is straightforward 
to show that the resulting plasmonic modes retain linear 
dispersions having associated velocities and Luttinger pa- 
rameters given by 



2 




_£± 

VF 



K,3 



vf/v± 



(38) 



where g is the coefficient of the {dxOc-){dx<j>c+) term and 
is given hy g = 4e^ft-3C/6/7rK. The associated tunneling 
density of states into the bulk of the tube is given by 



1 

Ctbulk = TT 



K, 



(39) 

This expression differs from the usual form (see the ex- 
pressions below Eq. ([M)) and below Eq. ([M]) ) because of 
the unusual charge-current coupling term in the interac- 
tions. 

Figure [5] shows the dependence of Kc- on electric field 
for a (20, 20) nanotube (with = 1/5) in the presence 
of a 5 T magnetic field. Interestingly, the fields render 
Kc- larger than unity, reflecting tendencies towards per- 
fect conduction in the c— sector. Furthermore, the values 
attained by Kc- show significant deviation from unity, 
making it feasible to observe the proposed spin-charge- 
band separation for these field configurations. Figure [S] 
shows a plot of Kc+ for the same situation. 



IV. LUTTINGER LIQUID PHASES 

The field-tuning of the Luttinger parameters discussed 
in the previous section offers a viable way of tuning the 
groundstatc of the nanotube through different phases and 
different ordering tendencies. In the absence of fields. 
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Egger and Gogolin [tII performed an involved analysis us- 
ing renormalization group arguments, refermionization 
and considerations of various susceptbilities to predict 
the ordering tendencies of the nanotube as a function 
of Luttinger parameters and temperature. In particu- 
lar, in the 'Luttinger liquid regime' which is the eas- 
iest to experimentally access, wherein all four sectors 
c± and s±, remain ungapped, the prediction is that 
for the range of interaction values Kc+ > 1/5, the sys- 
tem tends to show an inter-sublattice spin-density wave 
(SDWtt) ordering while for Kc+ < 1/5 it tends to show 
inter-sublattice charge-density wave (CDWtt) ordering, 
where the corresponding operators are defined, respec- 
tively, as OcDWl ^ J2paa i'laai' 

and OsDW 

J2paa (^'^laa'^-p±a<y; ■ Thcsc aualyscs involved consider- 
ing operators associated with certain orderings and de- 
termining the slowest decaying, equivalently, the most 
relevant operator (i.e., of the smallest scaling dimension). 

Here, we discuss the key changes that occur in the Lut- 
tinger regime in the presence of fields. We focus on the 
Luttinger liquid regime and consider the manner in which 
the field-tuned change in Luttinger parameters affect var- 
ious susceptibilities. We only consider the cases where 
either only an electric or magnetic field is present; the 
cases when both fields are present are extremely involved 
and beyond the scope of this paper. We do not take into 
account the effect of non-quadratic bosonic terms gen- 
erated by the fields; even if relevant, we expect that the 
bare coupling associated with these terms is so small that 
they only come into play at very low temperatures and 
not in the Luttinger liquid regime. 

Case of B = 0, E 7^ - For the case of only an electric 
field present, as discussed above, the effect of the field 
goes purely into changing the value of Kc+ ■ Given that 
experimentally the value of Kc+ is around and oftentimes 
higher than 1/5 and that the field tends to reduce the 
value of A'c+, the electric field provides a unique means 
of tuning from a tendency towards (SDWtt) ordering to 
that of (CDWtt) ordering (see for example figure [6]). 

Case ofB^^O, E = - The case of only a mag- 
netic field present, as discussed above, presents a slightly 
more complex situation in which both Kc+ and Kc- 
deviate from unity. As a result, various susceptibil- 
ities acquire a Kc- dependance in their scaling be- 
havior. For instance, operators associated with intra- 
sublattice ordering such as OcDwn ~ Y.paa i'laai' 

and OsDWo ~ J2paa '^'^laai'p-acr, which in the absence 
of fields are marginal, both acquire a scaling dimension 
(A'c- + K^_} + 2)/4. Tendencies for superconducting or- 
der become slightly stronger, though still irrelevant; the 
singlet pairing operator Osca ~ Y^paa crippaaipp-a-a ac- 
quires the scaling dimension (A'c- + K^_^ + 2)/4. 

To determine which ordering dominates, we consider 
the most relevant candidates: the CDWtt and SDWtt 
operators, both of which have scaling dimension {Kc- + 
Kc++2) / A; parts of the second order CDWtt operator de- 
noted by O^j^yyi that have scaling dimension Kc-+Kc+; 



and a fourth order CDWtt operator denoted by Ocdwtt 
which has scaling dimensions 4:Kc+. Comparing these 
exponents shows that Ocdwtt and Osdwti are more rel- 
evant than Oqj^^y^ for 15Kc+ > 2 + Kc-, a condition 
easier to satisfy in the presence of fields than in the field- 
free case since Kc- can then be less than 1 . Now OQ]^^/^ 
is more relevant than Ocniv-n- ^^^r 3Kc+ > Kc~. For 
^cdwtt to more relevant than OcDWir and Osdwh 
requires Kc+ + Kc- < 2/3, a condition requiring inacces- 
sibly strong interactions. Finally, to determine whether 
Ocdwtt or Ospw-K dominates, we appeal to the argu- 
ments of Ref. [3] ; at lower temperatures where the physics 
is dominated by certain strong coupling fixed points, pin- 
ning of the 9a+ mode suppresses Ocdwtt, making its 
magnitude in the Luttinger phase smaller than that of 
Osdwti- Though the methods of refermionization em- 
ployed to reach this conclusion are no longer valid for 
arbitrary values of Kc- ^ 1 , the strong coupling analysis 
still holds and we believe that a similar conclusion can 
be reached for the finite magnetic field situation. 

To summarize our results, while a detailed analysis, 
and considerations of operators that are not taken into 
account in Ref. p] are in order, the most likely scenario is 
that the nanotube in the Luttinger regime for the B ^ 0, 
E = case is dominated by SDWtt ordering tendencies 
for 15Arc-i- > 2 -t- Kc- (the more likely scenario gives that 
the deviation of Kc- from unity is not very large) and 
CDWtt ordering tendencies for 15A'c+ < 2 + Kc-- 



V. QUANTUM DOT PHYSICS 

For high resistance contacts or sufficiently low temper- 
atures [3| , the nanotube shown in Fig. [1] is only we akly 
coupled to the leads, thus forming a quantum dot [2l|. 
The behavior of such dots and related Coulomb block- 
ade effects have been extensively studied by theory and 
experiment 0, |23|. Typical of quantum dot physics. 
Coulomb blockade peaks have provided information on 
single-particle level spacings and charging energies associ- 
ated with the dot, as well as phase shifts due to scattering 
at the edges of the nanotube dot. Moreover, under cer- 
tain conditions, the nanotube quantum dot has revealed a 
periodicity of four associated with the degeneracy emerg- 
ing from the band and spin degrees of freedom. Recent 
work has also investigated the effects of a transverse mag- 
netic field on quantum dot behavior and the associated 
single-particle and charging energies [Tsj . Here, we study 
the role of boundary conditions on nanotube quantum 
dot physics, which requires subtle considerations in the 
presence of fields. We then discuss the quantum dot be- 
havior described by a finite-sized version of the nanotube 
Luttinger liquid Hamiltonian which takes into account 
relevant boundary conditions. 
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A. Field dependent single-particle energy spectrum 

As is well known, finitcncss of the tube length leads 
to a quantized single particle spectrum which depends 
on the boundary conditions at the ends of the tube. We 
assume that the wavefunctions at a given end are related 



by i'Racr = J2c 



where M is a matrix 



which depends on the the microscopic details of the tube 
end but is assumed to be energy independent. We spe- 
cialize to the case that the boundary conditions do not 
affect spin; that is we take Maaaa' = Saa'Saa', where 
5aa' is the Kronecker delta function. We thus assume 
the absence of magnetic impurities and local moments at 
the tube ends. 

In order to obtain the appropriate boundary condi- 
tions for the case of an asymmetric dispersion, we de- 
mand that the first quantized kinetic energy operator 
Hq = —iti'^j-a fVrdx together with the boundary condi- 
tions is self-adjoint [1^. This treatment does not account 
for the effect of interactions on the boundary conditions 
which would be more naturally discussed in terms of the 
bosonic fields. Such an analysis shows that there is an 
additional term in the current proportional to 52 ~ .94 
which vanishes for the density-density interaction con- 
sidered here [13] ■ By definition. 



Ho*) = -ifiYl J rVri^lo^da^Aa- (40) 



Since the boundary effects are assumed to be independent 
of spin, we have dropped the spin index. For an arbitrary 
spinor ^ with ^ = {tpR iPl)'^ where V'ii and are both 
two-component spinors in the Fermi point basis {ipr = 
{ipr+ i'r-)^), self-adjointness gives 

(*,i7o*) = (i?o*,*). (41) 

Integrating the left-hand side of this equation by parts 
gives 



cx=± 

Q = ± 



Self-adjointness is satisfied as long as the boundary terms 
vanish, and this leads to 



(42) 



with y^VR/vL S unitary. 

The details of the S-matrix can vary for each exper- 
imental set-up and depend on physical attributes such 
as the substrate, the hardness of the confining potential 
offered by the leads and the orientation of the tube's sub- 
lattices with respect to the leads. These parameters can 



be incorporated as variables in the boundary conditions 
which can then be utilized to obtain the single-particle 
spectrum. The most general version of these boundary 
conditions are outlined in Ref. [2^ via the effective-mass 
model. 

For a given S-matrix the spectrum of single particle 
states can be determined by applying the condition of 
Eq. at both ends and demanding that both the left 
and right movers have the same energy. The two Fermi 
points give rise to two sets of bands. The energy between 
two adjacent states in the same band is equal to TrhvH/L 
where vh is the harmonic mean, 



VH 



2vrVl 
vr + vl' 



(43) 



However, the energy offset of the bands from the Fermi 
energy depends on the details of the S-matrix. In general, 
the two Fermi points will give rise to two sets of energy 
states given by ■nnhvu I L -f Ai and nnfivH /L + A2 where 
neZ. 

Since we are ultimately interested in the spacing be- 
tween Coulomb blockade peaks, we focus on the energy 
difference between bands which we define as A^and = 
Ai — A2 (and for convenience we define Aband such that 
|Af,and| < TrhvH/L). Wc examine two special cases for 
the S-matrix; deriving the Aband for the most general 
scattering matrix would be a straightforward extension. 
First, consider the case in which the tube ends do not 
mix the Fermi points, though we allow the phase shift the 
electron suffers at the tube end S± (x) to be different for 
the two Fermi points (a = ±) and the two tube ends {x = 
0,L). In this case we have 5'++(a;) = \/vL/vRe^^+^^\ 

S—{x) = ^vlIvr^^-^'^'i and 5-+(a;) = S+-[x) = 0. 
The energy offset between the bands takes the form 



^band — -T 



2tA& 



27r Httvh/L 



(44) 



where 5i = {S+{L) - 6+{0)) - {6^{L) - 6^{0)) and 
J-"(x) = x— [xj and [x\ is the greatest integer less than or 
equal to x, and the quantity tAs is the field induced off- 
set between the two Fermi points as defined in Eq. (|10p . 
(See Appendix [C] for derivation.) 

Now, consider an electron that is completely scattered 
into the opposite Fermi point at both boundaries. For 
simplicity we take Saa'{x) = \/vL/vRfi^^'^^^ 5a,-a' ■ In 
this case, the splitting between bands takes the form 



T:hv 



H 



^band 



L 



-kpL 



2tAs 
nhvH I L 



f VR 


- vl\ 


\VR 





(45) 

For the limiting case of no fields (this also means that 
Vh ~ vr = vl = vp), one expects the existence of sets 
of four single-particle states, namely two degenerate sets 
of spin states and two sets of band states whose energy 
splitting depends on the various phase shifts and the ex- 
tent to which modes at the two Fermi points mix. For no 
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Fermi point mixing, the interband splitting is hvpSi/2L 
while for complete Fermi point mixing, the splitting is 
jF[2£fcF/7r]. Coulomb blockade experiments have 
shown an interband band splitting of about 10% [l3| of 
the Trhvp/L. Such a persistent approximate degeneracy 
in band energies for a range of tubes [l^ suggests that 
the magnitude of Fermi point mixing in these samples is 
minimal. 



B. Tunability of energy subband splitting 

As discussed above, the boundary conditions in a given 
experiment are not directly observable since Ahand de- 
pends on several parameters. Fields provide a way of 
controlling Aband as well as studying its physical origin 
in a particular sample. By scanning through various field 
strengths, the variation of the band offset can reveal in- 
formation about the nature of boundary scattering. For 
example, the extent to which a given tube interpolates 
between the two expressions given in Eqs. and ([^S]) 
can be used to determine the importance of (Fermi point) 
backscattering at the tube ends. In the case of a natu- 
ral band degeneracy in a tube (that is, no Fermi point 
mixing at the ends of the tube and 6i =0), both elec- 
tric and magnetic fields need to be applied to break the 
degeneracy; the magnitude of the subband splitting as 
a function of fields can be extracted from Eq. by 
setting Si ~ 0. An alternative approach for breaking the 
subband degeneracy was explored by Ref. [1^ in which 
a nonuniform external potential along the tube was ap- 
plied. However, this approach becomes infeasible for the 
case of a diagonal scattering matrix since it relies on band 
curvature away from half-filling. 

Thus combining electric and magnetic fields can pro- 
vide a means of breaking and tuning the degeneracy of 
the quantum states of electrons inhabiting the nanotube 
quantum dot. Of the four possible states discussed above, 
where the direction of spin is defined with respect to 
the magnetic field, an extra electron would occupy the 
ground state, which can be chosen to be any of the four 
depending on the direction of the fields. The quantum 
state can be characterized by a superspin inhabiting a 
SU{2) (g) SU (2) band and spin space. The enhanced con- 
trol of the spectrum of nanotubes that fields offer would 
obviously have important implications for any potential 
quantum information applications. 

C. Coulomb blockade physics 

We now consider quantum dot phenomena by taking 
into account interaction effects in addition to the single- 
particle level spacing analysis of the previous subsections. 
To present a coherent picture, we work within the context 
of the Luttinger liquid description for field configurations 
that retain gapless modes. This approach neglects the 
exchange energy within a dot, which while shown to be 



present, is often much smaller in magnitude than the level 
spacing and interaction energy . While our treatment 
captures salient features of quantum dot behavior, a full 
analysis of the Luttinger liquid formulation for the most 
general boundary conditions is yet to be performed. 

Following the method of Ref. Q (see also Refs. [13, 
H^), for a finite sized version of Eq. ([501) . decom- 
pose the bosonic fields 9a and tpa into sums of topolog- 
ical modes corresponding to a net occupation num- 
ber of the 'a' sector Na = /g dxO^ dx and harmonic 
modes corresponding to plasmons. For simplicity, we 
consider the case of no Fermi point scattering so that the 
Fermi point basis as defined by Eqs. ((2T|) and ([24]) and 
the band basis which diagonalizcs the boundary condi- 
tions in the previous subsection coincide (and therefore 
take a = c/s±). Assuming the boundary conditions de- 
rived in the previous section, we integrate out the x- 
dependence for the topological sector in the finite sized 
version of Eq. ((30)) . The resulting Hamiltonian associated 
with 'charging energy' for each topological sector takes 
the form 

-Ha = leaN^ (46) 

where = 4- AEa) and Ea is equal to the interac- 

tion energy of a given mode. The interaction energy for 
the net charge sector comes from the forward scattering 
contribution of Eq. ([25]) to yield Ec+ ~ V{k)/L. The 
contribution due to the electric field given by Eq. ([77]) 
is found to be fourth order in fields and can thus be ne- 
glected. The magnetic field, however, does contribute to 
the charging energy; from Eq. ((^51) , it can be shown that 
Ec- = 2e^h2b'^ / kL. This expression represents an upper 
bound that assumes the limit of no Fermi-point mixing. 

The topological modes correspond to the addition of 
charges onto the dot; in addition, plasmon modes that 
correspond to harmonic vibrations of the densities in the 
various sectors arepresent. In principle, the procedure 
employed in Ref. p, [23| to derive the structure of these 
plasmonic excitations can be generalized to the case of 
the asymmetric dispersion by incorporating the asym- 
metric description in Ref. [l9|. Here we forego such a 
derivation; most quantum dot experiments involving adi- 
abatic tuning of parameters such as gate voltage and thus 
probe purely ground state properties determined by the 
topological sectors. 

In the presence of a gate voltage Vq, the Hamiltonian 
associated with the topological modes of the nanotube 
quantum dot is given by 

Hl= J2 '^<'- ^^'=+ + \^bandN,. - AzN,+ , (47) 
a—±c/ s 

where ^ is essentially eVc and the term Az = hbB ac- 
counts for the Zeeman splitting. The hierarchy in energy 
scales can be summarized for a typical tube length of 
L = 500 nm. The intrasubband splitting is 3.3 meV. 
Therefore we have that ec_ w ^s+/- = 3.3 meV. The 
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charging energy is then ec+ = 47.7 meV. Consider a tube 
for which the ends of the tube do not appreciably change 
the Fermi point. In that case the sign of A(,and can be 
changed by reasonable electric and magnetic field values 
(for example, this is true for an n = 15 tube, a 6 T mag- 
netic field and an electric field of order 1 V/nm). Finally, 
fiB ~ 0.058 meV/T and therefore for most situations the 
Zeeman splitting will be at least an order of magnitude 
smaller than the other effects considered here. 

The Hamiltonian in Eq. (j47p forms a starting point for 
analyzing quantum dot and Coulomb blockade behavior 
in short nanotubcs. Typically, conductance is measured 
across the tube as a function of the applied gate volt- 
age Vg and a bias voltage Vb; while for the most part, 
energetic costs impede the fiow of electrons, at special 
degeneracy points that equally favor an occupation of N 
and N + 1 electrons, zero bias Coulomb blockade con- 
ductance peaks can be observed. Given that the occupa- 
tion numbers Na with a = {c/s, ±} are good quantum 
numbers, the net energy of the system = {Hl) for 
a given configuration of electrons corresponds to a given 
combination of eigenvalues of Na- The equilibrium con- 
figuration of electrons on the quantum dot can thus be 
derived by minimization, i.e., by requiring = for 
all Na sectors, subject to the physical constraint that 
that electron occupation numbers Na=±,a-='j/i take on 
integer values. The relationship between these two bases 
is given by Nc± = (A^+t + N+i) ± (iV_| + N^i) and 
Ns± = {N+^ - iV+i) ± (iV_T - iV-i). 

As an illustrative example, consider the Coulomb 
blockade situation for the simple case of no fields, no 
band or Zeeman splitting (Ahand = A^ = and only the 
charging energy in the c-\- sector is non-zero). Applying 
the condition that |^ = gives ec+Nc+ — 4^ = and 
Nc- = Ns+ = Ns- = (where fj, is varied by Vg). These 
conditions suggest that all the states at a given energy 
level will fill before filling the subsequent energy level. 
Now, as the chemical potential is increased, the first extra 
electron added onto the dot can occupy any of the degen- 
erate states of the \a — ±,(j =t / j) space. Suppose that 
this electron goes in the a = -I- band with its spin up. 
The state of the tube is then characterized by the quan- 
tum numbers {Nc+, N^-, Ns+, N^^) = (1,1,1,1). Fur- 
ther increasing the chemical potential by an amount Ec+ 
will add the next electron to any of the remaining three 
states. For example if the filling obeys Hund's rule (i.e., 
assuming the exchange interaction which we have thus far 
neglected) then band a ~ — would be filled with spin up 
electron thus giving the state (2,0,2,0). An additional 
third electron can occupy -|- i or — | state. The energy 
cost for adding each of these extra electrons is Ec+, re- 
fiecting the Coulomb charging energy. However, the fifth 
extra electron requires a chemical potential increase of 
+ Ec+, reflecting the Coulomb energy as well as 
the excitation energy required to occupy the next energy 
level. This analysis, executed within the Luttingcr liquid 
description, re pro duces the periodicity of four observed 
in experiment [l0[. 



In the presence of fields and intrinsic subband splitting, 
an analysis similar to the one above can be performed for 
altered values of e and KaS and the orders of magnitude 
discussed after Eq. (|i7)) . We take Az « and A^and > 
0. Furthermore, suppose that for a given tube that a 
nonzero magnetic field gave Ec- < Af,and/2. In that 
case, the tendency to minimize Nc- would give rise to a 
shell-filling opposite to the usual Hund-likc filling. For 
example, the order in which states are filled could take 
te forn - t, - i, + + i or equivalently (1, -1, 1,-1) 
(2, -2, 0, 0) ^ (3, -1, 1, 1) ^ (4, 0, 0, 0). Instead of the 
Coulomb blockade peaks being of periodicity 4 described 
above (spaced apart by Ec+, Ec+, Ec+j Ec+ + ttHiih/L, 
spacings become Ec+ + Ec- , Ec+ - iE^- + A^and, Ec+ + 

Ec-,Ec+ + Ec- + TThvH/L - Aband- 

The effect of shell-filling in nanotube quantum dots has 
been investigated experimentally by Liang et al. [l3| . The 
four-electron periodicity was clearly observed via trans- 
port measurement. Parameters of charging energy and 
exchange energy were determined. The above results, in- 
cluding field-dependent Luttinger parameters, the effect 
of boundary conditions, the band splitting, asymmetric 
dispersions and additional two-electron periodicity can 
thus be studied in a similar setup with transverse fields. 
Additionally, the Fabry-Perot transmission resonances in 
the presence of a transverse magnetic field predicted by 
Bellucci and Onorato [l^ , as well as possible resonances 
by both transverse electric and magnetic fields (and their 
relative angles), can be investigated in the same experi- 
mental setup. While our arguments here have been con- 
fined to adiabatic tuning and zero-bias conductance, our 
approach can be used to investigate non-equilibrium phe- 
nomena, temperature dependences, higher order tunnel- 
ing events such as cotunneling and non-adiabatic tuning. 
Each of these considerations, which is beyond the scope 
of this paper, would involve excitations of the plasmonic 
modes. 



VI. DISCUSSION; RELEVANCE TO 
EXPERIMENT 

We have investigated the effects of transverse electric 
and magnetic fields on armchair carbon nanotubes. We 
found that fields can break several symmetries inherent 
to the carbon nanotubes — the valley degeneracy, the left- 
right-mover degeneracy, and the particle-hole symmetry. 
The magnitude of a gap in the nanotube spectrum can 
be continuously tuned by varying the strength and the 
relative angle of the two fields. We also found that the 
electron-electron interaction is modified by both fields 
and thus Luttinger-liquid parameters can be tuned by 
fields. In particular, an interesting consequence is the 
possibility of spin-charge-band separation. We also dis- 
cussed how the fields can be used to study boundary 
effects in finite sized tubes and to describe the Coulomb- 
blockade physics in presence of fields. 

Each of these salient features can become manifest in 
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experiment, some in dramatic ways. While we summa- 
rize these experimental signatures here, details of the 
physics and orders of magnitude can be found in the 
relevant section. At the band structure level, the reduc- 
tion of the Fermi velocity can be observed by measuring 
the particle level spacing in a finite size tube. The shift 
in the Fermi momentum induced by either field may be 
detected by virtue of the associated Friedel oscillations 
around a dopant or impurity using a STM (see for ex- 
ample [2§|) For both fields present at an arbitrary angle 
to one another, a continuous conduction gap occurs at 
the Fermi energy, discernible via direct conductance mea- 
surements, shifts in Coulomb blockade peaks and STM 
measurements [3, BIIH- In the transport measurement, 
the conductance peak can vary from to 2e^//i, and to 
4e^//i, depending on the fields and the chemical poten- 
tial. Perhaps the most dramatic prediction of band struc- 
ture analysis is that an electric and a magnetic field per- 
pendicular to one another and the tube axis would give 
rise to a current carrying ground state. Thus, a tube 
subject to this field configuration placed across two leads 
should induce a measurable current even in the absence 
of an applied voltage drop across the leads. 

In the regime in which SWNTs exhibit Luttinger liquid 
behavior, strong enough fields can give rise to significant 
changes. The value of the Luttinger parameter associated 
with the net charge sector can be tuned via either an elec- 
tric field or a magnetic field or both. Furthermore, the 
presence of an electric field gives rise to density-density 
interactions associated with the difference in densities in 
the two bands and results in the deviation of the asso- 
ciated Luttinger parameter Kc- from its non-interacting 
value of unity. The tunneling density-of-states, a quan- 
tity ubiquitous to a range of experiments including scan- 
ning tunneling microscope studies and conductance mea- 
surements 0, m HO, im , would reflect these changes in its 
power-law behavior. 

As discussed in Scc lIVl changing the Luttinger pa- 
rameter via fields can result in tuning through phases 
having charge-density-wave order or spin-density-wave 
order; such phases are potentially measurable in STM 
and neutron scattering experiments. Another exciting 
prospect comes about in the tuning of Kc- , namely, that 
of spin-band-chargc separation. In the past, momen- 
tum resolved tunneling experiments have resolved charge 
and spin modes moving at different velocities in quantum 



wires [20| : in principle, a more elaborate version of such 
an experiment could detect charge, spin and band modes 
moving at three different velocities in nanotubes. 

In the quantum dot regime, the application of fields 
acts as a controlled means of changing the Coulomb 
blockade structure of the dot, and could potentially have 
a plethora of applications. For tubes that preserve the 
four-fold degeneracy emerging from spin and band de- 
grees of freedom in the absence of fields, the presence of 
fields can serve to break this degeneracy. For tubes that 
show a lack of degeneracy, fields provide a way of deter- 
mining the origin of degeneracy breaking. This effect has 
potential applications to quantum information. Through 
achieving a desired amount of degeneracy breaking for 
the four states, fields can be an effective means of ini- 
tializing the quantum state of an extra electron added 
on to the dot. Having initialized a quantum state, tran- 
sitions can be induced to other states. Additionally, as 
has been demonstrated for semi-conducting quantum dot 
spin states [s^l, superpositions can be created by apply- 
ing oscillating fields. For the energies quoted above, os- 
cillation frequencies would be on the order of 10^^ Hz. 

To conclude, transverse fields induce a rich range of 
physical effects in the electronic properties of SWNTs 
from band structure effects to one-dimensional behavior 
to quantum dot physics. Several of these features are 
very much within the reach of current experimental ca- 
pabilities and are of both fundamental and applied value. 
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symmetries of the band structure and Jeremy McMinis 
for his assistance in preparing the nanotube figure (Fig. 
1). This work was supported by the NSF under the grant 
DMR-0605813 (SV,WD). TCW was supported by IQC, 
NSERC and ORF. 



APPENDIX A: BAND STRUCTURE 
CALCULATION 

In the absence of any fields, the eigenstates of an in- 
finitely long armchair nanotube are superpositions of the 
states l^A/B^) defined by Eq. (1). For the particular 
case of an armchair nanotube Eq. (4) gives 



^il-^sl^s) = -— ^e^2.^^<jeW +2e"^cos 



2 h \2tt J 



2n 27r 
cos —s — cos — 



2^/3 



where h = "^fj^l^^ ■ For smaU magnetic fields (6 <C 1) for £ - £' = ±1 mod N and 
we have . 



{^\\Hb\'^'b) =tb sin 



m-e- ( 



(Al) 



(A2) 
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Applying perturbation theory near the Fermi points 
requires care because of the near degeneracy of the left 
and right movers. In this case, the usual procedure of 
first diagonalizing the nearly degenerate subspace fails 
because the matrix elements within this subspace vanish 
to the order we are working. However, there are matrix 
elements for transitions to other energy levels {£ ^ 0) 
and these matrix elements give rise to an effective inter- 
action between states in the nearly degenerate subspace. 
The subspace can be diagonalized once these additional 
interactions are taken into account [1^. 

For the gapless case, the low energy spectrum near 
half-filling is 



e^„(fc) = hTVF,r \ k - akFj + atAs + 0{k'^), (A3) 
The renormalized Fermi velocity is given by 

Vr = vf{1- Avib^ - Av2Uy ± AvsbUy) , (A4a) 



where 



Avi 



AV2 



5 -I- 4 cos - 

I n 

3(l + 2cos^)^ 



t- cos 4- 2 cos ^ 

3n 3n 

12 (1 



(A4b) 



- cos - 

n 



i) 



(A4c) 



AV3 



cos 



6n 



cos 1^) CSC TT- 

on / on 



^^(l + 2cos^)^ 



and 
kp 



3a 



2V3(l + 2cos^) V3 



sm 



6n 



(A5) 

The shift between the two Fermi points (see Eq. (10)) is 
given by 



%/3sin 
As = — ^rWy. 

1 + 2 cos :f ^ 

3n 



(A6) 



For the case of mutually perpendicular fields discussed in 
section II, the electronic densities are given by the vector 

1 / I + g^ru cos {% + -^]-g2ba cos + ^ 
^™^^ 2 1 - .gi^cos 92bacos - ^ 

(A7) 

where the upper and lower components are the electronic 
densities over the A and B sublattices respectively. The 
constants gi and 52 are given by 



1 TT 

32 = VScsc— . 

2n 



(A8) 
(A9) 



APPENDIX B: INTERACTION TERMS 

In order to find the form of the effective interaction 
^aa' need account for not only the radial dependence 
of the wave functions but also for the physical separa- 
tion between the sublattices. Although the factorization 
performed in Eq.[T2]is an approximation, we may still ac- 
count for the physical separation of the sublattices. We 
follow the approach introduced in Q . 



2nR r27rR 



ds ds' 



JO (27ri?)2 



where U{d) is a shorthand for the Coulomb interaction 
with an offset d, that is 



U{d) = U{x-x',s^s' + d) 



(Bl) 



where the right hand side of this equation is given by Eq 
The constant cq which appears in Eq. (j26p is given by 



Co(n) = -7 - ^J^ rf^ln[cos - + — ] . 

Similarly we find that the values of hi and /12 defined 
in section 3 are given by 

(A4d) hi{n) = {cl - cl) h{n) + 2ciC2/2(n) + 2 (c? + d) h{n), 



h2{n) = {cj - 4) h{n) + 2c3C4/2(n) + 2 (c^ + cl) h{n), 
and 

/i3(n) = 2[(ciC3 -I- C2C4) fi{n) + [cia + C2C3) f2{n) 

-|-2(ciC3 + C2C4)/3(n)], 



where 
h{n) 



ttB. 



dz 
7vR 27ri? 



In 



2R \2RJ 



In 



^ 2R \2r) 



cos|. 



f2in) 



ttB. 



dz 
ttR 27ri? 



In 



Qc / y 



2R \2RJ 



In 



2 z + ac / y 

' 2R '^\2r) 



h{n) 



TiR 



dz 
2^ 



■In 



2R \2R) 
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and 



This gives constraints on the four momenta 



Cl 
C2 

C4 



1 TT TT 

- CSC — COS — 

2 2n 6n 

1 TT TT 

- CSC — sin — 

2 2n 6n 



= V3 



CSC 



TT 


TT 


TT 






COS — 




on 


on 


TT 


TT 


TT 




sm — 


sm — 


2^ 


6n 


3n 



The values of cq, hi, h2 and /i3 have been tabulated for 
various tube sizes and are displayed below. 



n 


Co 


hi 


h2 


h3 


5 


-0.239 


0.491 


0.112 


3.780 


10 


-0.064 


0.697 


0.040 


10.370 


15 


0.005 


0.786 


0.020 


17.417 


20 


0.025 


0.834 


0.122 


24.605 


25 


0.434 


0.865 


0.008 


31.855 



{kR±-kL±)L = 27rn± + 5± [L) - S± (0) , (C3) 



where n± are arbitrary integers. For the a = + branch, 
the energy levels corresponding to kfi^ and k^^ should 
be equal (in order for the wavefunction to represent an 
eigenstate). Hence 



hvR{kR+ - kp) + tAs = -hvL{kL+ - kp) + tAs, (C4) 
which gives 



kL+ = -kp - vr 



2mi+ + 6+iL) -S+{0) 
{vr + vl)L 



(C5) 



and thus 



APPENDIX C: SINGLE-PARTICLE SPECTRUM 
WITH BOUNDARY SCATTERING 

Here we illustrate the case in which the tube ends do 
not mix Fermi points; that is S+-^-{x) = \/vL/vRe^^+'^^\ 
S—{x) = y/vlj^e'^-^''^ and S^+{x) = S+-{x) = 0. 
Assume the wavefunctions arc 



27m+ + M£)-M0) ^ A (ra^ 
e+{n+)^hvH — hAs, (C6) 



where vh = "^vrVl/ {vr + vl)- Similarly, for a = — , we 
have 



(Cl) 
(C2) 



The energy spectrum for each branch is described by 
Eq. Uni), e.g., eR.± = hvR{kR± =F kp) ± tAs. The self- 
adjointness condition (|^^ at a: = 0, L gives 




e'-^-W J\d) 

gi(5+ (L) Q 

e'^-^^'> 



, , ^ 27rn_-F(5_(L)-(5_(0) . 
t_(n_)^hvH ^-As. (C7) 



Hence, we arrive at the interband energy difference 



Aband = e+(n+) -e_(ri-) = — - — [n- 



L \ 2tt hirvH / LJ ' 

where Si = [5+{L) - 6+{0)^ - (6-{L) - 5_(0)) and the 

h± are selected so that |A(,a„(i| is less than the intraband 
spacing. Hence wc obtain Eq. (HH). A similar considera- 
tion leads to Eq. ((45)) . 



[1] S. lijima. Nature 354, 56 (1991). 

[2] T. Giamarclii, Quantum Physics in One Dimension. Ox- 
ford Science Publications, (2004). 

[3] H. Ajiki and T. Ando, J. Phys. Soc. Japan 65, 505 (1996). 

[4] S. Saito, M. Dresselhaus, and G. Dresselhaus, Physical 
Properties of Carbon Nanotubes (Imperial College Press, 
1998). 

[5] S. lijima and T. Ichihashi, Nature 363, 603 (1993); D. 



S. Bethune, C. H. Klang*, M. S. de Vries, G. Gorman, 
R. Savoy, J. Vazquez, and R. Beyers, Nature 363, 605 
(1993). ' 

[6] C. Kane, L. Balents, and M. P. A. Fisher. Phys. Rev. 

Lett. 79, 5086 (1999). 
[7] R. Egger and A. O. Gogolin, Eur. Phys. J. B. 3, 281 

(1998). 

[8] M. Bockrath, D. H. Cobden, J. Lu, A. G. Rinzler, R. 



15 



E. Smalley, L. Balents, P. L. McEuen, Nature 397, 598 
(1999). 

[9] Z. Yao, H. W. C. Postma, L. Balents, and C. Dekker, 

Nature 402, 273 (1999). 
[10] W. Liang, M. Bockrath, and, H. Park, Phys. Rev. Lett. 

88, 126801 (2002). 
[11] U. C. Coskun, T.-C. Wei, S. Vishveshwara, P. M. Gold- 

bart, and A. Bezryadin, Science 304, 1132 (2004). 
[12] S. Sapmaz, P. Jarillo-Herrero, L. P. Kouwenhoven, and, 

H. S. J. van der Zant, Semicond. Sci. Technol. 21, S52 

(2006). 

[13] H. Ajiki, T. Ando, J. Phys. Soc. Japan 62, 1255 (1993). 
[14] D. S. Novikov and L. S.Levitov, e-print 

|arXiv:cond-mat/0204499 vl. 
[15] H.-W. Lee and D. S. Novikov, Phys. Rev. B. 68, 155402 

(2003). 

[16] S. Bellucci and P. Onorato, Eur. Phys. J. B 47, 385 

(2005) . S. Bellucci and P. Onorato, Ann. of Phys. 321, 
934 (2006). 

[17] N. Mason. Private Communication, (2008). 

[18] S. Bellucci and P. Onorato, Eur. Phys. J. B 52, 469 

(2006) . 

[19] M. Trushin and A. L. Chudnovskiy, EPL 82 17008 

(2007) . 

[20] O. Auslaender, A. Yacoby, R. de Picciotto, K. W. Bald- 
win, L. N. Pfeiffer, and K. W. West, Science 295, 825 
(2002). O. M. Auslaender, H. Steinberg, A. Yacoby, Y. 
Tserkovnyak, B. I. Halperin, K. W. Baldwin, L. N. Pfeif- 
fer, and K. W. West, Science 308, 88 (2005). 

[21] M. A. Kastner, Artificial Atoms, Physics Today 46, 24 



(1993). 

[22] Y. Oreg, K. Byczuk, and B. Halperin, Phys. Rev. Lett. 
85, 365 (2000). 

[23] Thanks to M. Stone for a helpful discussion of the bound- 
ary conditions and the condition of self-adjointness in this 
context. 

[24] V. I. Fernandez, Ph.D. thesis, Universidad Nacional de 
La Plata, 2002. 

[25] E. McCann and V. Fal'ko. J. Phys. Cond. Matt. 16, 2371 
(2004). 

[26] A. Makarovski, L. An, J. Liu, and, G. Finkelstein, Phys. 

Rev. B 74, 155431 (2006). 
[27] M. Fabrizio and A. O. Gogolin, Phys. Rev. B 51, 17827 

(1995). 

[28] S. Eggert, H. Johannesson and A. Mattsson. Phys. Rev. 

Lett. 76, 1505 (1996). 
[29] P.T. Sprunger, L. Petersen, E.W. Plummer, E. Laegs- 

gaard, F. Besenbacher, Science 275 1764 (1997). 
[30] M. Bockrath, D.H. Cobden, J. Lu, A. G. Rinzler, R.E. 

Smalley, L. Balents, and P.L. McEuen, Nature 402, 273 

(1999). 

[31] H. W. C. Postma, M. de Jonge, Z. Yao, and C. Dekker, 

Phys. Rev. B 62, RIO 653 (2000). 
[32] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 

K. C. Nowack, T. Meunierl, L. P. Kouwenhoven and L. 

M. K. Vandersypen, Nature 442, 766 (2006). 
[33] L. D. Landau and E. M. Lifshitz, Quantum Mechanics. 

Butterworth-Heinemann, (1958). 



16 



